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Introduction 



Existing proposals for next-generation Earth-based gravitational wave detectors, including advanced interferome- 
ters (IFOs) and wideband (nested) acoustic detectors, are aimed at constructing large-scale prototypes, with strain 
sensitivity goals of the order of 10~ 23 in a wide (a few kHz) band. The successful achievement of these goals will 
depend on substantial advances to be made, well beyond the limits of presently achieved figures of merit, in several 
critical areas, including (ultra)cryogenics (mK cryostats with adequate refrigerating power), material science (achiev- 
ing extremely high mechanical quality factors at cryogenic temperatures) and technology (annealing/sinthcrizing 
high-quality deca-ton size objects), electronics (SQL amplifiers), lasers and optics (high-finesse, large-waist, short 
optical cavities performing as displacement readouts at sensitivity levels ~ 10 -23 m Hz -1 / 2 ). All these goals are two 
order of magnitude beyond the best figures available today, on average. 

The MAGO proposal stems from a different perspective, and emphasizes a relatively new concept. The MAGO 
design aims at constructing a relatively large (N > 10) number of detectors, relying only on available and already 
proven technology. 

Experience gained on downsized MAGO prototypes suggests that MAGOs will be relatively cheap, compact-sized, 
and structurally simple. Indeed, within the limits of available and already proven technologies, the proposed MAGO 
design could allow to construct compact detectors, featuring strain sensitivities comparable to those of present day 
cryogenic acoustical detectors, at a fraction of their size and cost. This would permit to place several MAGOs even in 
a small area (which would be required, e.g., for observing a stochastic GW background), either co-tuned or frequency 
staggered, thanks to the unique bandwidth and center frequency tuning ease offered by the design. 

The potential advantages of many-detector arrays and networks have been repeatedly emphasized. Similar to 
radioastronomy, gravitational wave astronomy will be eventually made possible by the availability of large assemblies 
of co-tuned (array) or frequency staggered (xylophone) detectors featuring robust, high duty-cycle operation, as a 
result of constructive simplicity, and non critical core technologies. The MAGO project and design is definitely aimed 
at moving a first, perhaps tiny step toward this direction. 

The following specific technical points are worth being emphasized: 

• The individual MAGO bandwidth and center frequency of operation can be easily tuned well beyond the 
operational limits of present day detectors, imposed by laser shot noise (IFOs) and Brownian (acoustic) noise; 

• The very large parametric-conversion gain which provides the basic principle of operation of the detector brings 
the equivalent temperature of the front-end electronics down to the standard quantum limit already using 
standard, cheap and reliable off-the-shelf HEMT amplifier (T ~ 2 K at 1 GHz), without the need of resorting 
to more sophisticated devices; 

• The required electrical quality factor for the Nb-coated superconducting cavities, has been safely and routinely 
achieved in accelerator technology, and could be realistically improved by one order of magnitude relatively 
soon; 



• The required mechanical quality factor of the cavities is also conservative and could be improved; 

• The phase noise of the microwave pumping signal is well within the limits of currently available technology. 

The main technology challenge is thus essentially related to the design and testing of a cryostat capable of removing 
the heat produced in the cavity walls by the microwave pumping signal (a few Watts, typically) so as to maintain an 
operating temperature ~ 2 K, while guaranteeing adequate acoustical isolation and very low intrinsic noise. 

On the other hand, it should be noted that the MAGO design has no a-priori intrinsic limitation which could 
prevent it from reaching performance figures, both in terms of strain sensitivity and bandwidth, comparable to those 
of proposed next generation large-scale detectors (including advanced IFOs and wideband nested acoustic detectors) , 
provided one pushes the relevant critical design figures (cryogenics, mass, mechanical quality factor, readout noise) 
up to comparable levels. 

In particular, if we were able to design and operate the cryostat down to the mK range, the MAGO strain sensitivity 
could easily reach the h ~ 10~ 23 level, while preserving MAGO's almost unique features in terms of center frequency 
and bandwidth easy tunability. As of today, the best performing lOOmK cryostats are capable of delivering only a 
few hundred mW. We note that next generation optically read-out wideband nested-acoustical detector would face 
the same technological challenge, among others, in order to reach their foreseen sensitivity goal. 

I. EXPERIMENT OVERVIEW 

In the last decades, several laboratories all around the world have promoted an intense effort devoted to the direct 
detection of gravitational waves. The detectors, both those in operation and those being developed, belong to two 
conceptually different families, massive elastic solids (cylinders or spheres) Q and Michelson interferometers 0. Both 
types of detectors are based on the mechanical coupling between the gravitational wave and a test mass, and in both 
types the electromagnetic field is used as motion transducer. 

In a series of papers, since 1978, it has been studied how the energy transfer induced by the gravitational wave 
between two levels of an electromagnetic resonator, whose frequencies ui\ and L02 are both much larger than the 
characteristic angular frequency Q of the g.w., could be used to detect gravitational waves 03- The energy transfer 
is maximum when the resonance condition \u>2 — Wi| = fi is satisfied. This is an example of a frequency converter, i.e. 
a nonlinear device in which energy is transferred from a reference frequency to a different frequency by an external 
pump signal. 

In the scheme suggested by Bernard et al. the two levels are obtained by coupling two identical high frequency 
cavities Each resonant mode of the individual cavity is then split in two modes of the coupled resonator with 
different spatial field distribution. In the following we shall call them the symmetric and the antisymmetric mode. 
The frequency difference of the two modes (the detection frequency) is determined by the coupling, and can be tuned 
by a careful resonator design. An important feature of this device is that the detection frequency does not depend 
on its mechanical properties (dimensions, weight and mechanical modes resonant frequencies), though, of course, 
the detector can be tuned so that the mode splitting equals the frequency of a mechanical resonant mode. The 
sensitivity in this and other experimental situations will be discussed in the following. Since the detector sensitivity 
is proportional to the electromagnetic quality factor, Q, of the resonator, superconducting cavities should be used for 
maximum sensitivity. 

An R&D effort, started in 1998, was completed at the end of 2003 00. Its main objective was the development 
of a tunable detector of small harmonic displacements based on two coupled superconducting cavities. Several cavity 
prototypes (both in copper and in niobium) were built and tested, and finally a design based on two spherical cells 
was chosen and realized (Fig. The detection frequency, i.e. the frequency difference between the symmetric and 
antisymmetric modes, was chosen to be: u>2 — ~ 10 kHz (the frequency of the modes being Wig » 2 GHz). An 
electromagnetic quality factor Q > 10 10 was measured on a prototype with fixed coupling. 

The tuning system was also carefully studied. The coupling strength, and thus the tuning range, is determined by 
the diameter of the coupling tube and by the distance between the two sperical cells. A central elliptical cell, which 
can easily be streched and squeezed, was found to provide a tuning range of several kHz (4-20 kHz in the final design) . 
A prototype with the central elliptical cell was built and tested (Fig. 

The system was also mechanically characterized, and the mechanical resonant modes in the frequency range of 
interest were identified. In particular the quadrupolar mode of the sphere was found to be at 4 kHz, in good 
agreement with finite elements calculations. 

The detection electronics was designed. Its main task is to provide the rejection of the symmetric mode component 
at the detection frequency. A rejection better than 150 dB was obtained in the final system. 

Starting from the results obtained in the last six years, we are now planning to design and realize an experiment for 
the detection of gravitational waves in the 4-10 kHz frequency range. Our main task is the design and construction 





FIG. 1: Artistic view of the coupled spherical cavities with the central tuning cell 




of the refrigerator and of the cryostat (including the suspension system), which houses the coupled cavities. The 
refrigerator must provide the cryogenic power needed to keep the superconductiong cavities at T » 1.8 K (approx. 10 
Watts) without introducing an excess noise from the external environment. A design based on the use of subcooled 
superfluid helium is being invesigated. 

In the following a detailed description of the various issues aforementioned will be given. Expected system sensitivity 
will also be discussed. 



II. PHYSICS MOTIVATION 



The spectrum of gravitational waves of cosmic origin targeted by currently operating or planned detectors spans 
roughly 1 from 10~ 4 to 10 4 Hz. 

The / < 10 _1 Hz region of the GW spectrum, including galactic binaries Q, (super)massive BH binary inspirals 
and mergers compact object inspirals and captures by massive BHs 0, will be thoroughly explored by LISA [To| . 
which might be hopefully flown by year 2015. Ground based interferometers and acoustic detectors (bars and spheres) 
will likewise co-operate in exploring the / > 10 1 Hz region of the spectrum, including compact binary inspirals and 
mergers |Tlj , supernovae and newborn black- hole ringings |l2j| , fast-spinning non-axisymmetric neutron stars jl.'il| , and 
stochastic GW background (l4|. 



1 We leave out deliberately the ELF (10 — 18 — 10 -15 Hz) radiation resulting from inflation-enhanced primordial gravitational fluctuations, 
expected to show up in the polarization anisotropy of the cosmic microwave background (M. Kamionkowski, A. Jaffe, Int, J. Mod. Phys. 
A16, ff6, 200f), the VLF (10 -7 — 10 — 9 Hz) radiation possibly resulting from extremely massive black-hole systems and early— universe 
processes (A.N. Lommen and D.C. Backer, Astrophys. J., 562, 297, 2001, astro- ph/0107470 I, exotic electromagnetic-to-gravitational 
wave conversion mechanisms, which might originate gravitational radiation in the VHF to SHF bands (Fang-Yu Li and Meng-Xi Tang, 
Int. J. Mod. Phys Dll, 1049, 2002) and the relic gravitational radiation (B. Allen and R. Brustein, Phys. Rev. D55, 3260, 1997). 



The whole spectral range from 1CP 4 — 10 4 Hz, however, is far from being covered with uniform sensitivity. Plans 
are being made for small-scale LISA-like space experiments (e.g., DECIGO, aimed at covering the frequency 
gap 10 _1 — 10 1 Hz between LISA and terrestrial detectors. 

Several cryogenic/ ultracryogenic acoustic (bar) detectors are also operational, including ALLEGRO , AURIGA 
[13, EXPLORER O, NAUTILUS 19], and NIOBE |oj. They are tuned at ~ 10 3 Hz, with bandwidths of a few 
tens of Hz, and minimal noise PSDs of the order of 10 -21 Hz -1 / 2 . 

Intrinsic factors exist which limit the performance of both IFOs and acoustic detectors in the upper frequency 
decade (/ > 10 3 Hz) of the spectrum. 

The high frequency performance of laser interferometers is limited by the oc f 2 raise of the laser shot-noise floor. 
While it is possible to operate IFOs in a resonant (dual) light-recycled mode, for narrow-band increased-sensitivity 
operation, the pitch frequency should be kept below the suspension violin- modes |2lj |. typically clustering near and 
above ~ 5 • 10 2 Hz. 

Increasing the resonant frequency of acoustic detectors (bars, spheres and TIGAs), on the other hand, requires 
decreasing their mass M. The high frequency performance of bars and spheres is accordingly limited by the cx M ~ x / 2 
dependance of the acoustic detectors' noise PSD. 

The next generation of resonant detectors will be probably spheres or TIGAs (Truncated Icosohedral Gravitational 
Antennas, |22j) 2 The MINIGRAIL |2jj spherical prototype 3 experiment under construction at Leiden University 
(NL), as well as its twins planned by the Rome group |24| and at Sao Paulo, Brazil |2^|, is a relatively small (CuAl 
(6%) alloy, 65cm, 1.15ton) spherical ultracryogenic (20mK) detector with a 230Hz bandwidth centered at 3250Hz, 
and a (quantum limited) strain sensitivity of h ~ 4 • 10 -21 . Spherical (or TIGA) detector might achieve comparable 
sensitivities up to / ~ 4 • 10 3 Hz. 

Summing up, the GW spectrum below / ~ 10 3 Hz might be adequately covered by ground-based and space-borne 
interferometers. The range between 10 3 Hz and ~ 4 • 10 3 Hz could be sparsely covered by new-generation acoustic 
detectors. The high frequency part (/ > 4 • 10 3 Hz) of the gravitational wave spectrum of cosmic origin is as yet 
completely uncovered. Within this band, GW sources might well exist and be observed 4 . Indeed, the ultimate goal of 
gravitational-wave astronomy is the discovery of new physics. In this spirit, the very existence of gravitational wave 
sources of as yet unknown kind could not be excluded a-priori. 

The above brings strong conceptual and practical motivations for the MAGO proposal. The MAGO design is 
easily scalable, and may be constructed to work at any chosen frequency in the range 10 3 — 10 4 Hz, with uniform 
(narrowband) performance. On the other hand, the MAGO instrument appears to be comparatively cheap and 
lightweigth, thus allowing to build as many detectors as needed to ensure adequate covering of the high frequency 
(/ ~ 4- 10 3 Hz) GW spectrum. In view of their limited cost, MAGOs might also be nice candidates for many-detector 
networks, to achieve very low false alarm probabilities in coincidence operation. In addition MAGO-like detectors 
operating at / ~ 10 3 Hz might hopefully provide coincident observations with both acoustic detectors and IFOs, 
based on a different working principle. 

Before all this might come into reality, it will be necessary to build and operate one or more MAGO prototypes so 
that some basic issues might be efficiently addressed and solved, viz.: 

• efficient decoupling from platform — > suspensions design; 

• efficient and quiet cooling to 1.8K — > cryostat design; 

• efficient readout — > microwave feeding and tapping networks, and low noise amplifier design. 

In parallel, a start-to-end simulation codes should be implemented, in order to tune all design parameters for best 
operation. In particular, criteria for obtaining the best tradeoff between detector bandwidth and noise levels should 
be investigated, with specific reference to selected classes of sought signals. 



2 Spheres and TIGAs share the nice feature of being inherently omnidirectional, and should allow to reconstruct the direction of arrival 
and polarization state of any detected gravitational wave, by suitably combining the outputs of transducers gauging the amplitudes of 
the five degenerate quadrupole sphere modes (C. Zhou and P.F. Michelson, Phys. Rev. D51, 2517, 1995). 

s A cryogenic solid-CuAl sphere resonating at ~ 10 3 Hz would be ~ 4 m and weigh in excess of 100 tons. The related technological 
challenges could be alleviated by resorting to hollow geometries (J. A. Lobo, Class. Quantum Grav., 19, 2029, 2002). 

4 A well known back-of-an envelope estimate (motion at the speed of light along body-horizon circumference) gives the following upper limit 

for the spectral content of gravitational waves originated in a process involving an accelerated mass ~ M: f sup < 4 C GM ~ 10 4 (Mq/M) 
[Hz]. 



III. SOURCES IN THE RANGE 10 3 HZ-10 4 HZ 



A. Mergers 

A well known upper frequency limit for the spectrum of gravitational radiation originated in a process involving an 
accelerated mass M is given by 

corresponding to the rather extreme assumption of motion at the speed-of-light at the gravitational body horizon. 
These conditions are met almost verbatim in compact binary mergers, which have been accordingly indicated as 
the most reliable sources of gravitational waves in the frequency range from 10 3 Hz to 10 4 Hz [54(. Gravitational 
waves from mergers are expected to carry rich astrophysical information on the nature of the coalescing stars, and in 
particular, as to whether exotic (e.g. quark) matter is involved 

In view of our still relatively limited knowledge of the relevant waveforms > h is difficult to estimate the related 
signal to noise ratios. However mergers will follow an inspiral phase, for which more reliable estimates can be made. 

To learn as much as possible from the waves generated during the merger, broadband GW detectors should be 
supplemented by narrow band detectors. A "xylophone" of narrow band detectors would probe features of the merger 
waveform which should be robust in the sense that they would not require detailed modeling of the waves' phasing. 
To make best use of these detectors, the network should be designed in an optimal way: the narrowband detectors 
should be tuned, in concert with the broadband detectors, so that the network of all detectors is most likely to provide 
new information about merger waves Practical considerations such as cost, available facility space and tunability 
strongly favor the use of a network of MAGOs. 



1, Event rates 



Crude estimates of NS-NS merger event-rates in our Galaxy are obtained by dividing the estimated number of 
NS-binaries in our Galaxy by the average time they take to merge. Multiplying this figure by the number of galaxies 
(or fraction of Galaxy disk volume) within the reach (visibility distance) of a given detector gives a crude estimate of 
the detectable event-rate. 

This exercise has been repeated through the last decade by several Authors |H^-|Hflj. 

The most refined analysis available to date, including the simulation of selection effects inherent in all relevant 
radio pulsar surveys, and a Bayesian statistical analysis for the probability distribution of the relevant merger rate, 
has been presented by Kim, Kalogera and Lorimer |^|. They estimate the NS-NS merger event-rate in the Galaxy 
between 0.2 x 10 -5 Mpc~ 3 yr -1 and 6.0 x 10 -5 Mpc~ 3 yr _1 depending on the assumed pulsar population model, and 
provide a NS-NS merger event-rate probability distribution yielding a most likely value for the NS-NS merger rate 
between 0.3 x 10 -5 Mpc~ 3 yr _1 and 1.7 x 10 -5 Mpc -3 yr _1 events per year at a 68% confidence level, in agreement 
with previously obtained orders of magnitude. 

Interest in mergers on behalf of the GW Detectors' Community has been revived by the recent (published dec. 
2003) discovery by Burgay and co-workers of a new, highly relativistic NS-NS binary (PSR J0737-3039) with an 
orbital period of 2.4 hr and an estimated lifetime of only 85 Myr [6(J. The estimated NS-NS merger event-rates 
depend critically on the shortest observed NS-binary lifetime. This observation alone suffices to boost the estimated 
NS-NS merger event rate in the Galaxy by a factor between 4 and 7 at the 68% confidence level, as shown in Fig. 

m 

The remarkable weakness of PSR J0737-3039 radio signals, despite its relative nearness (~ 600 pc) to Earth, further 
suggests that there might be a plethora of as yet undetected short-lived NS-NS binaries in the Galaxy. This is the 
dominant source of uncertainty in estimating NS-NS merger rates, which might accordingly be ~ 10 2 larger than 
estimated, as stressed in |6lj| . where it is argued that all other modeling uncertainties build up only to an uncertainty 
factor of the order of unity. 

In order to translate the above event rates into observable event rates, one has to estimate the visibility distance of 
the MAGO detector. We define here the visibility distance as the distance of a source yielding SNR= 1. 

To determine the visibility distance of MAGO (and/or a MAGO array, or xilophone), one may adopt the std. 
(conservative) procedure appropriate for the case where detailed knowledge of the waveform is not available. One 
accordingly has to compare the squared sought signal characteristic amplitude, 

2 _ 2 dE 

c ~ ^~dj 
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FIG. 3: Probability density function (pdf) for the increase in the double-neutron-star merger rate resulting from the discovery 
of PSR J0737-3039. The dashed line represents the pdf corresponding to the old merger rate (before the discovery of PSR 
J0737-3039). The heavy solid line represents the pdf corresponding to the new merger rate (from 



to the squared detector's r.m.s. noise ampltude, 

h 



2 _ fS(f) 



(F 2 

so as to recast the signal to noise ratio into the form 



SNR: 



h 2 (f) ll/2 
d(log/)- c '" 



K(f) 



where S(f) is the detector noise power spectral density (PSD), F is the detector's directivity function, (■) denotes 
averaging over all directions, and dE/df describes the spectral energy content of the sought signal. Based on reasonable 
assumptions about the transition from inspiral to merger, we might use for this latter the formula (in geometrized 
units) 

dE ^ f 0.91 M 2 X (4/x/M) 2 , / G (f merge, fring) 

df \ 0, elsewhere 

first proposed in where f me rge ~ /Isco ~ 4 x 10 3 (A/ Q /M) Hz and f Hng - 2.6 x 1O 4 (M /M) Hz is the 
quasi-normal mode ring-down frequency of the final collapsed object |6^ . 

For a medium-size MAGO {Sl /2 - 6 x 10- 21 Hz" 1 / 2 @ 2 KHz in a 350 Hz bandwidth) the corresponding visibility 
distance (SNR= 1) will be ~ 0.2 Mpc for a M - 10Af Q merger with a sharp cut-off at M ~ 2M Q . With a MAGO 
xylophone covering the 2-8 KHz band at the same PSD level, the mass range could be extended down to M ~ 0.6M Q . 
(see fig. EJ). 
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FIG. 4: Visibility distance of one MAGO (dashed line) and of a xylophone of six MAGOs covering the 2-8 kHz frequency range 
(solid line) vs. total mass of the merger (in solar mass units). 



B. Bursts 



The characteristic noise amplitude of one MAGO is h n (f c ) ~ 7 x 10 19 at the characteristic frequency f c ~ 2 kHz. 
This accounts for the isotropic conversion into GW energy of 8 x 10~ 4 solar masses located at a distance of 8 kpc 
(SNR=1) (see fig. EJ. 

Recently [H3 a coincidence excess was found among the data of the resonant bars EXPLORER and NAUTILUS, 
when the detectors are favorably oriented with respect to the galactic disk. The observed coincidence corresponds 
to a conventional burst with amplitude h c ~ 2 x 1CP 18 and to the isotropic conversion into GW of 4 x 10~ 3 M Q , 
with sources located in the Galactic Center. Still unknown phenomena other than GWs, though, cannot be ruled out 
as causes of the observed events A possible way to distinguish GWs from other sources is the use of several, 

different, detectors with good sensitivity. From this point of view the use of a network of MAGOs could contribute 
in clarifying this subject. 



C. MACHO Binaries 



The gravitational wave frequency corresponding to the last stable circular orbit, loosely marking the transition from 
inspiral to plunge in binary coalescence, is roughly |fi3 |: 

~ 3976^ [Hz] (3.2) 

for a (symmetric, non spinning, circular orbit) binary with total mass M. 

The spectral window of MAGOs would thus be highly appropriate to observe gravitational waves from black-hole 
MACHO binaries jjj^j with a typical total mass of 0.6M Q < M < 2M Q . It is speculated |H3 that one half of the 
galactic halo mass (corresponding to ~ 3 x 10 11 M Q ) consists of MACHOs with a typical mass of M ~ Q.3M®. Among 
these, ~ 30% should be bound in binary systems, formed at the same time (10 10 yr) the galaxy was born and uniformly 




FIG. 5: Visibility distance of one MAGO vs. burst energy (in solar mass units). 



distributed across the halo. The rate of MACHO binary coalescence could be of the order of 10 events per year 
in the galactic halo. 

A nice representation of current expectations for BH-MACHO binary inspiral (and merger) observable event rates 
has been recently given by J.C.N, de Araujo and co-workers 68] for the (advanced) Brazilian spherical antenna "Mario 
Schenberg" as a function of the burst sensitivity, and is shown in fig. The burst sensitivities h in fig. |H|are readily 
translated into visibility distances d by noting that h{d\) / h(d,2) — dij&\. 

For visibility distances up to the Galaxy's border (d ~ 20 Kpc) the event rate increases linearly up to ~ 0.05 yr . 
This accounts for the linear segments on the right side of fig. From the border of the Galaxy up to the outskirts 
of M31 and M32, at d ~ 700 Kpc, intergalactic BH-MACHOs might be visible. Their contribution, however, would 
be negligible under the assumption that they follow the distribution of dark matter. This corresponds to the plateau 
in figure [fl] As the visibility distance overpasses Mil and M32, the event-rate contribution of these latter is added, 
corresponding to the step in fig. H3 For larger distances, beyond 1 Mpc, the background BH-MACHOs contribution to 
the event rate becomes eventually dominant. The relevant event-rate is obtained using the coalescency rate computed 
in (leftmost region in fig. HJ. 

The visibility distance (SNR=1) for observing a coalescing O.5M© + O.5M MACHO binary with a MAGO [S^ 2 = 
6 x 10~ 21 between 1850 Hz and 2150 Hz) would be ~ 18 kpc. The maximum visibility distance for this detector is 25 
kpc for observing a 0.9M Q + 0.9M Q MACHO binary 5 . 

Furthermore, in view of MAGO's relatively simple, lightweight and hopefully cheap features, one might think of 
building arrays of many MAGOs, either identical or differently tuned detectors (xylophone) whose response would 
allow to deduce the relevant source parameters (chirp-mass etc.). Fig. shows the visibility distance of both six 
identical MAGOs as a function of the binary total mass. 



5 The visibility distance of the initial "Mario Schenberg" (aka mini-GRAIL, SFERA) antenna (S* /2 ~ 2 X 10~ 21 Hz" 1 / 2 @ 3.2 kHz in a 
50 Hz bandwidth) will be d ~ 17 kpc at SNR=1 for a O.5M + 0.5M Q BH-MACHO. 




h 

FIG. 6: BH-MACHO binary inspiral (and merger) observable event rates for the (advanced) Brazilian spherical antenna "Mario 
Schenberg" as a function of the burst sensitivity (from [6£| )• 



D. Submillisecond Pulsars 



The existence of sub-millisecond pulsars has been speculated by various Authors. It is well known that the (lower) 
limiting spin period of a neutron star depends on the assumed equation of state, being ~ 600 fis for the softest one [jjj . 
Submillisecond pulsars, if any, could be strange (quark) stars |7l| as well. Known millisecond pulsars are believed to be 
"recycled" neutron stars spun-up by accretion from the companion, in low-mass X-ray binar y sy stems (LMXB) [72j . 
On the other hand, the so called gravitational radiation induced Rossby(r)-mode instabilities j73| are now believed to 
play a crucial role in setting an upper limit to the spin period 74]. As a matter of fact, the distribution of periods 
of all known pulsars drops off quit e sharply at about 2 msec. All sub-millisecond pulsar search campaigns made so 
far have been indeed negative |75| . |76j |. Under such circumstances the possible direct observation of gravitational 
radiation from pulsars above 2 KHz appears to be extremely difficult, though it cannot be completely ruled out. 



E. Stochastic Background 

The standard procedure to search for a stochastic background of gravitational waves (unresolved superpositions of 
gravitational wave signals of astrophysical or cosmological origin) is to cross-correlate (and suitably filter) the data of 
two (or more) detectors 6 

Specific strategies for GW background detection using two bar detectors [8(| , two spheres [8lJ , two interferometers 
|82j . a bar and an interferometer [83| have been discussed by various Authors, and are in current use to set upper 



The power of a statistical hypothesis testing strategy based on a single detect or's output to decide whether any excess noise could be 
interpreted as due to such a stochastic background would be exceedingly poor l77l . 
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FIG. 7: Visibility distance of an array of six identical MAGOs vs. the binary total mass. The typical distance of the galactic 
halo is also shown (gray area). 



limits on the gravitational wave background power spectral density on the basis of available data from operational 
gravitational wave antennas |84|-|85|. 

We plan to develop a similar (straightforward) analysis for the special relevant cases of two MAGOs, one MAGO 
and one bar (or sphere), and one MAGO and one interferometer. 

In view of the easy tunability of MAGOs at any frequency in the range between 10 3 Hz and 10 4 Hz and beyond, 
using MAGOs might be interesting for probing the otherwise unaccessible high frequency part of the sought stochastic 
gravitational wave spectrum. 



IV. DETECTOR LAYOUT 
A. Electromagnetic design 

In order to build an efficient detector, a suitable cavity shape has to be chosen. According to some general 
considerations, a detector based on two coupled spherical cavites looks very promising (Fig. The choice of 

the spherical geometry is based on several factors. From the point of view of the electromagnetic design the spherical 
cell has the highest geometric factor G, thus it has the highest electromagnetic quality factor Q, for a given surface 
resistance R s (Q = G/R s ). For the TEon mode of a sphere, the geometric factor has a value G w 850 0, while 
for standard elliptical radio-frequency cavities used in particle accelerators, the TMoio mode has a value G w 250 0. 
Looking at the best reported values of surface resistance of superconducting accelerating cavities, which typically are 
in the 10~ 8 f2 range, we can extrapolate that the electromagnetic quality factor of the TEon mode of a spherical 
superconducting cavity can be Q s» 10 10 — 10 11 . 

In the first generation of detectors, dedicated to the development of the experimental technique, the internal radius 
of the spherical cavity will be r *=s 100 mm, corresponding to a frequency of the TEon mode u> ~ 2 GHz. The overall 



system mass and length will be M « 5 kg (with a wall thickness w ~ 2 mm) and 1/ w 0.5 m. The choice of the wall 
thickness is made considering both practical and design constraints. On one hand, the wall thickness should be kept 
small enough to allow an easy fabrication while maintaining sufficient stiffness to withstand the external pressure 
once the cavity is evacuated. Furthermore, wall thickness was chosen to optimize the cavity cooling process, and to 
guarantee optimum stability against point-like thermal dissipation due to possible defects present on the cavity inner 
surface. On the other hand, wall thickness can be used to design a particular mechanical resonant frequency and it 
is obviously related to the mass of the detector, which plays an important role in the signal to noise ratio. 

Since this type of detector is ideally suited to explore the high frequency region of the g.w. spectrum, we plan to 
build a tunable cavity with 4 kHz< u>2 — u>i < 10 kHz, which is outside the spectral region covered by the resonant and 
interferometric detectors, both existing and planned, and is still in a frequency region where interesting dynamical 
mechanisms producing g.w. emission are predicted 27, 28, 29] . 

The interaction between the g.w. and the detector is characterized by a transfer of energy and angular momentum. 
Since the helicity of the g.w. (the angular momentum along the direction of propagation) is 2, the g.w. can induce a 
transition between the two levels provided their angular momenta differ by 2; this can be achieved by putting the two 
cavities at right angle or by a suitable polarization of the electromagnetic field inside the resonator. The spherical 
cells can be easily deformed in order to induce the field polarization suitable for g.w. detection. The optimal field 
spatial distribution has the field axes in the two cavities which are orthogonal to each other (Fig. |SJ . Different spatial 
distributions of the e.m. field (e.g. with the field axes along the resonators' axes) have a smaller effect or no effect at 
all. 




FIG. 8: Electric field magnitude of the TEon mode. Note the alignment of the field axis 

A tuning cell is inserted in the coupling tube between the two cavities, allowing to tune the coupling strength (i.e. 
the detection frequency) around the design value. The dependence of the detection frequency on the distance between 
the two coupled cells is shown in Fig. while its dependence on the diameter of the tubes is shown in Fig. ITUl 

One detector based on two spherical niobium cavities (with fixed coupling) has recently been built and tested at 
CERN. A second detector with variable coupling has also been built and is now being tested (Fig. 

The first test on the cavity with fixed coupling, showed a quality factor Q > 10 10 (see Fig. ITT|) . This corresponds to 
a surface resistance R s ~ 50 nf2, a factor of ten higher than the best values reported for superconducting accelerating 
cavities. The obtained result is very satisfactory. In fact, the whole fabrication procedure (including surface treat- 
ments) is optimized for the elliptical cavity geometry used for high energy particle acceleration. Some development is 
still needed to tailor the technique to the spherical shape of our resonator and to obtain a surface quality comparable 
to that routinely obtained on elliptical cavities that would lead to a quality factor Q w 10 11 . 

B. Mechanical design 

From the mechanical point of view it is well known that a spherical shell has the highest interaction cross-section 
with a g.w. and that only the quadrupolar mechanical modes of the sphere do interact with a gravitational wave 30]. 
The mechanical design is highly simplified if a hollow spherical geometry is used. In this case the deformation of the 
sphere is given by the superposition of just one or two normal modes of vibration and thus can be easily modeled. 
In fact, the proposed detector acts essentially as an electro-mechanical transducer; the gravitational perturbation 
interacts with the mechanical structure of the resonator, deforming it. The e.m. field stored inside the resonator 
is affected by the time-varying boundary conditions and a small quantity of energy is transferred from an initially 
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excited e.m. mode to the initially empty one. We emphasize that our detector is sensitive to the polarization of the 
incoming gravitational signal. Once the e.m. axis has been chosen inside the resonator, a g.w with polarization axes 
along the direction of the field, will drive the energy transfer between the two modes of the cavity with maximum 
efficiency. With standard choices for the axes and polar coordinates, the pattern function of the detector is given by 
F x = — cos(0) sin(20), and is equal to the pattern function of one mechanical mode of a spherical resonator. 

C. Status and results on the simulation of the mechanical modes 

We have run an ANSYS mechanical simulation of our prototype over a wide frequency range (0-7000 Hz), finding 
approximately 180 resonant modes and we have compared these modes to the (analytical) modes of a spherical shell 
for identification. 

We foresee to suspend the detector by a region near its center of mass, a region where the tuning system should 
also be implemented. For this reason, the model used for the simulation consists of half detector (single cell) fixed 
on the central pipe, near the whole detector center of mass. Let us focus on two frequency ranges: the modes below 
the quadrupole frequencies (ANSYS calculated v qU ad = 3646 Hz, theoretical frequency is v qua d — 3750 Hz) and the 
modes above v qua d- 

Below Vq U ad, we find a very limited number of modes (~ 16), well separated in frequency. These modes are due 
to the movement of the pipes with respect to the cavity (bending, rotation and cantilever modes). Above v qua d we 
find a great number of modes, only related to the cell. Fig. 1121 shows the cumulative number of modes versus their 
frequency. The slope of the ANSYS calculation (blue dots) is very close to the slope of the analytical spherical modes 
(green cross), meaning that the number of the single cell modes found above v qua d corresponds to the ideal sphere, 
although several mode shapes may be distorted due to the presence of the pipes. 

The frequency range around the quadrupolar mode of interest (7V,L,M = 1,2,2) gives us two main hints: 

• The rightfulness of the detection mode, which is minimally perturbed by the pipes attached to the cavity and 
oscillates approximately at the wanted frequency; 



10 5 r-. — 1 — 1 — i — 1 — 1 — r 




1000 ' — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 

22 24 26 28 30 32 

Diameter [mm] 

FIG. 10: Detection frequency vs. coupling tubes diameter 



• The presence of only quadrupolar modes in the range of interest. 

The first quadrupolar mode resembles one of the fundamental quadrupole mode of a single, ideal sphere. The five 
quadrupolar modes are all near to each other, although the cavity pipes interference gives a discernible spread. A 
visual identification with the ideal spherical modes can be done for the L = 2, L = 3 and L = 4 mode groups. Higher 
order modes, although somewhat similar to the ideal case, are in general less easy to identify. 

1. Role of other mechanical mode scarce 

Consider now a detector tuned to an angular frequency f2 = u>2 ~ wi and let us focus on OMR operation, so that 
> ijJmech- The signal, driven by the gravitational wave, comes from the quadrupolar modes (in particular, for a 
well oriented source, from the N,L,M = 1,2,2 mode). That is, the gravitational wave couples to those modes, driving 
them at angular frequency f2. The thermal noise, which will surely come at least from the same quadrupoles, may in 
principle come from other mechanical modes. In particular, we seek deformations which contribute to the noise but 
are not affected by the signal (g.w.). For a mode to contribute to noise sources though, it must fulfill two stringent 
requirements: 

a) The coupling coefficient C^ 1 = J" (H2 • Hi — E2 • Ei)£ m • dS, between the electromagnetic modes and the 
mechanical mode must not be null; 

b) The mode frequency must be relatively near to the detection frequency. 
The first requirement is the most important. 

In the ideal, spherical case, the coupling coefficient C^ 1 is different from zero for the class of quadrupole modes alone. 
Other angular numbers exhibits a strictly zero coefficient due to symmetry mismatch between the electromagnetic 
field and the mechanical mode. 

A clue for this behavior comes from the chosen electromagnetic mode angular dependence, which is oc cos (29) (it 
is not depending on <fi), while the radial component of the displacement is proportional to cos (L9). This calculation 
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has been performed on an ideal, hollow sphere and using analytical expressions for the electromagnetic fields. We 
understand that these results must be checked against the more realistic deformations delivered by ANSYS, and our 
future efforts will surely tackle this problem. We are confident though that the general behavior will be consistent, 
since the preliminary analysis shows a considerable similarity between sphere modes and the ANSYS model. 

V. SUSPENSION SYSTEM 

Like a resonant bar, the whole detector must be thoroughly insulated from the environment and the mechanical 
design of the insulation is therefore strongly bonded with the cryogenic design. In particular, the insulation for the 
direct LHe supply to the cavity (see section IVIfl may prove to be a formidable task. The mechanical insulation 
foreseen for this detector may be inspired by the already existing resonant bars suspension systems, although some 
major differences exist. 

Since one of the feature of this detector is the tunability, the mechanical design should provide adequate insulation 
in the frequency range 4-20 KHz and a very rough estimation suggests a minimum value of 250 dB of attenuation in 
the desired band (by comparison, Al bars suspension are estimated to provide sa 300 dB of insulation). A sketch of 
the suspension system is shown in fig. 1131 

The room temperature filter bank is attached to the outer cryostat hull and will provide a first decoupling to 
the environment. There follows a cryogenic box, a facility which will host some critical low-temperature electronic 
components and the source for the superfluid LHe, and finally the cryogenic filter bank, which should provide the 
required insulation. This second stage filter is connected to the innermost cryostat, which is a superfluid LHe container 
hosting the detector and its tuning system. The non trivial task of filtering out environmental acoustical noise on 
a wide frequency band can be more easily accomplished by optimizing the two filter banks on separate frequency 
ranges, for example, the first stage can be optimized for low frequencies. The effectiveness of the cold filter (second 
stage) can be weakened by the superfluid LHe fccdthrough, which runs from the cryogenic box to the MAGO cavity 
cryostat. Although the superfluid LHe can be conveyed by means of a suitable capillary tube, acoustic vibrations in 
LHe, triggered in cryo-box for instance, require a deeper understanding and maybe a dedicated damping system. The 
mechanical insulation design is made more complex by the MAGO cavity plus the superfluid cryostat weight. Since 
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FIG. 12: Spectrum of the mechanical normal modes of the small-scale detector in the range < / < 7000 Hz. 



their weight is of the order of some tens of kilos, it is comparable to the suspensions weight and therefore, the load 
on uppermost suspension element can be considerably different from the lowermost one. All these problems will be 
addressed in the following years, where a comprehensive cryogenic and mechanical design will constitute the main 
effort. 



A. Seismic Noise 

Preliminary seismic noise measures have been performed in our laboratory during working hours (day time). The 
measures were taken in a range from 10 to 9000 Hz, which should cover the sensitive region for the MAGO detector. 
By comparing the data to the expected level of the thermal noise we can estimate the required attenuation of the 
external environment noise to be (at least) of the order of 180-200 dB in the range 500 Hz-9 kHz, and 160 dB if we 
limit the frequency range to / > 4 kHz. 

1. Data taking 

Seismic noise was acquired with a PCB 393B12 accelerometer connected to a Stanford research SR780 FFT ana- 
lyzer. The accelerometer has a nominal frequency range of 0.05 to 4000 Hz and a sensitivity of 1019.4 mV/(m/s 2 ). 
Although that particular accelerometer is not fit for measuring frequencies up to 9000 Hz, we have nevertheless taken 
the frequency range 10-9000 Hz in order to have a preliminary, order of magnitude measure of the environmental 
disturbances in the range of interest. We are justified in taking the high frequency measures because the accelerometer 
response has a rising characteristic with increasing frequency (ref. model 393B12 operating guide). It is thus expected 
that the high frequency portion of the seismic noise power spectral density might be overestimated. More accurate 
evaluation of the seismic noise in the high frequency range (/ > 4000 Hz) will be done as soon as our group is provided 
with a suitable accelerometer (for instance the PCB 307B). 

The accelerometer was secured to a 25 kg lead brick, which allowed us to easily sample the vibrations along the 
three orthogonal directions in space, defining axis z oriented perpendicular to the floor, whereas x and y parallel to 
the room walls. For each direction, a measure with the accelerometer switched off (null measure) was performed in 
order to take into account the SR780 FFT and cables intrinsic noise over the whole frequency range. 



First stage 
suspensions 



LHe, vacuum 
and signal 
feedthrough 





Cryo-box (superfluid 

LHe facility and 
critical electronics) 



Second stage 
suspensions 



Superfluid LHe 
cry o stat 



MA GO cavity 



Tuning system 



FIG. 13: Sketch of the suspension system 



The sub-cooled liquid helium environment, needed for MAGO operation, is provided by means of a big mechanical 
pump, located next to the laboratory. When switched on, this pump emits an audible noise, and for this reason 
measures were taken with the pump both switched on and off. 

The Volts measured as the output of the accelerometers were converted into the equivalent force acting on our 
detector. This allows direct comparison with the Langevin force spectral density and therefore yields immediately 
the required attenuation. Sticking to a conservative approach, we have assumed that the noises measured in the three 
directions equally couple to the detector. 



2. Results 



Figure IT"T1 shows the total force spectral density (summed over x, y and z) with the helium pump switched on and 
off. The contribution of the pump is relevant up to a frequency / ~ 3 kHz. 

If we require that the measured force spectral density be of the same order of the Langevin force spectral density 
Sff = AK])TMLj m /Q m , we immediately can calculate the approximate attenuation for the range of frequencies greater 
than a chosen cut off frequency / > /o (see fig. IT5|) . 

Although these measures are not conclusive, we expect the attenuation values proposed in this document to be a 
fairly good approximation. More accurate and reliable numbers could be provided with suitable test equipments and 
a setup nearer to a g.w. detector experiment. Note that the assumption of equal coupling between the vibrations in 
the three directions and the system may well be an overestimation of the real case. Further studies could therefore 
relax the attenuation requirements. 
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FIG. 14: Total force spectral density measured in the "quiet" condition (helium pump switched off, dotted line) and "noisy" 
condition (helium pump switched on, continuous line). 



VI. CRYOGENICS 

Assuming the conservative value of 1 W/m 2 for the super-insulation radiation losses, the radiating heat load is 
Prad ~ 4 W. Adding to this figure a couple of watts to account for the conduction losses in the cavity suspension 
(tie-rods), we obtain the value of Ptot ~ 6 W as an upper limit. These steady state losses have to be compared with 
the RF cryogenic losses of our detector. 

Assuming, in operation, a peak surface magnetic field H sur f = 0.1 T (half of the critical field of niobium) the stored 
energy in the cavity will be U « 30 J. We remark that H sur f = 0.1 T corresponds (in an accelerating cavity) to 
25 MV/m of accelerating field, a value routinely achieved in the cavities developed for the Tesla project. Given the 
geometric factor of our cavity, G — 835 f2, and assuming a residual surface resistance of 5 nQ (the surface resistance 
routinely obtained in the LEP-II cavities at 1.8 K), we can foresee a quality factor Q w 10 11 , corresponding to a 
dissipation of 4 W in the helium bath (@ H sur f — 0.1 T). 

The planned operating temperature of our system is 1.8 K to fully exploit the advantages of the RF supercon- 
ductivity. At the chosen frequency of 2 GHz the surface resistance at 1.8 K is well saturated at the residual value, 
avoiding the change of surface resistance produced by the heating of the surface. The use of superfluid helium as 
refrigerant guarantees a very good thermal dissipation for the cavity, with an even temperature distribution along the 
cavity surface. 

To further improve the heat exchange and reduce the effect of the helium boil-off we propose a scheme of refrigeration 
similar to the one foreseen for the LHC magnets and already used since the eighties for the refrigeration of high field 
magnets (e.g. TORE-II supra). The flow chart of the refrigerator is shown in figure [TBI 

This refrigeration scheme, using subcooled superfluid helium, at atmospheric pressure and 1.7-1.8 K, avoids all the 
problems related to the bubbles in the refrigerating bath, improves the Kapitza resistance at the helium-cavity wall 
interface and increases by 30% the peak transfer value to the helium bath. This refrigeration scheme will be the best 
solution also from the point of view of the mechanical noise induced on the detector by the LHe bath. The subcooled 
operation at 1.8 K and atmospheric pressure (compared with saturated helium at 1.8 K and 20 mbar), will eliminate 
the possibility of the bubble gas production. This effect of an abrupt reduction of the bath induced mechanical noise 
below the helium A-point was already experienced on Weber type gravitational detectors [3lJ . 
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FIG. 16: Flow chart of the subcooled superfluid helium refrigerator 



To form a bubble of a certain size in subcooled helium we need not only to transfer to the helium bath the amount 
of heat corresponding to the heat of vaporization, but also the enthalpy needed to reach the normal boiling point of 
the gas at 4.2 K. Furthermore the Claudet-type refrigerator reduces to a minimum the refrigerator-bath interaction 
minimizing the mechanical noise coming from the pumps and from the saturated Hel bath. The sub-atmospheric 
Hell bath cools the main bath only by conduction via the heat exchanger HT, the refill of this bath is done trough 
the JT needle valve working on the liquid helium flow. The hydraulic impedance of this needle is quite high and 
damps any pressure fluctuation due to the pumping system. The 4.2 K bath is in hydraulic contact with the 1.8 K 



sub-cooled bath via the high impedance duct to reduce the heat input to the superfluid helium and allows to reduce 
the bath temperature well below the A-point of the helium (2.19 K). To obtain this result we need a quite long 
and narrow hydraulic channel with a very high thermal impedance, but also the hydraulic impedance of the channel 
is high, decoupling the superfluid bath from the saturated (4.2 K) or nearly superfluid bath. In this way pressure 
fluctuations due to the He bath refilling or to turbulences, produced by thermo-acoustic oscillations of the saturated 
bath, will have a negligible effect on the superfluid bath. 

Coming to the conclusions the proposed refrigeration scheme will greatly help in reducing the sensitivity of the 
proposed detector to the acoustic noise produced by the fluctuation of the helium bath wile giving us the more 
comfortable situation by the point of view of the RF performances. 



A. Refrigeration scheme and Noise Issues 

The microwave power dissipation in the cavity walls is easily found by the quality factor definition 

Q=f (6.D 

with u> angular frequency of the radio-frequency (RF) field, U Electromagnetic energy stored in the resonator, and 
P RF power lost in the cavity wall. 

In our case (operating frequency u> ~ 2 GHz) for a stored energy U ~ 18 J, corresponding to cavity operation at 
peak surface magnetic field B pea k ~ 100 mT (already achieved on the prototype cavity measurements), the power 
dissipation is P ~ 20 W at the design electromagnetic quality factor Q = 10 10 . The measured quality factor for the 
prototype MAGO cavity is Q — 5 x I0 10 at the design field, giving a safety factor of five over the power dissipation 
level quoted in the proposal. 

The RF power dissipation is fairly constant on more than 70% of the resonator surface; moreover the temperature 
dependence of the RF surface resistance of the niobium is highly non linear [8^|, the thermal conductivity of the 
superconductor is very poor, and the heat capacitance low |8ft chapter 15]. For those reasons (for the quoted quality 
factor and field) the resonator can only be cooled via an helium bath. 

Localized thermal connections to an heat sink at low temperature, as the one used for cooling bars, lead to thermal 
runaway of the resonator over the critical temperature of the niobium, up to 20 to 30 K. Considerations on the thermal 
noise of the detector, and the temperature dependence of the RF surface resistance of superconductors suggest to 
operate the detector at the lowest possible temperature. Practical considerations on readily available refrigerators 
with cooling power exceeding 10 watt, restrict our choice to a refrigerator scheme using superfluid helium operating 
around T = 1.5-1.8 K. 

At the operating temperature of T = 1.8 K the largest practical superfluid helium refrigerators, designed for the 
operation of superconducting cavities, have a cooling power in the multi kW range and sup erfluid helium refrigerators 
developed for superconducting magnet operation have ~ 10 kW of refrigerating power 90]. The world largest dilution 
refrigerator, the Compass refrigerator at CERN, has a maximum cooling power of 1 W at 0.3 K and 10 mW at 80 mK 
[9ll |. Finally, considerations of sec. IVI B 21 on the impact of the coolant fluid on the mechanical quality factor of the 
resonator, strongly suggest to use as coolant superfluid helium at T ~ 1.6-1.5 K. 



1. Refrigerator Scheme 

The proposed refrigeration system is a pressurized superfluid helium bath operating at T — 1.5—1.7 K and P = 
10 5 Pa, based on the original Claudet system developed at CEA for the TORE II Supra superconducting Tokamak 
[92| . A conceptual design of the refrigeration system is shown in fig. IT7I 

The cavity is immersed in a pool of stagnant superfluid helium and is cooled by conduction (not convection) by the 
surrounding pressurized superfluid helium bath at T = 1.7 K and P = 10 5 Pa. 

The operation of the cryostat is the following: the reservoir (housing the cavity) is filled with liquid superfluid 
helium at one bar. The superfluid helium in the reservoir is kept at T = 1.7 K by the Hell-Hell heat exchanger 
cooled at 1.5 K by the saturated superfluid helium produced by pumping on the saturated Hell pot. The pressure 
on the reservoir is kept at one bar by the hydraulic connection (a narrow channel) with the saturated liquid helium 
(Hcl) at one bar and 4.2 K. Inside this narrow channel lies the Hel-Hell interface with a temperature drop of 2.5 K, 
and no pressure drop because no liquid helium flow takes place in the channel. 

The heat produced by the RF power dissipation on the cavity walls is transferred to the Hell-Hell heat exchanger by 
the large thermal conductivity of the superfluid helium. This large thermal conductivity, (together with the peculiar 
properties of the superfluid helium) keeps homogeneous the temperature of the refrigerant in the reservoir, preventing 




FIG. 17: Pressurized Superfluid Helium (Claudet type) refrigerator. 



any convection in the liquid and the generation of any flow-induced noise in the reservoir. Moreover at the rated 
dissipation of 20 W, the power flow per unit area is roughly 10 mW(cm)" 2 , a factor of 100 lower than the value for 
the onset of bubble formation in the helium bath at the cavity surface. Roughly speaking, to produce bubbles in 
pressurized Hell one needs to transfer to the bath the extra energy needed to heat the liquid from 1.7 K to 4.2 K. 

Both effects, tightly connected to the use of pressurized superfluid helium, will guarantee us against noise generation 
in the process of cavity refrigeration. 

2. Induced noise in the cavity reservoir 

The proposed refrigeration scheme is quite safe in keeping to a minimum the noise generated from the thermal 
dissipation in the liquid helium bath. An additional source of acoustic noise is the pressure fluctuation in the helium 
reservoir. 

Two main sources inducing acoustic waves in the helium are the pressure fluctuations on the saturated 4.2 K helium 
bath and pressure fluctuations induced by the pumping system on the saturated 1.5 K Hell pot feeding the Hell-Hell 
heat exchanger. 

The pump-related pressure fluctuations on the Hell pot are decoupled from the pressurized Hell reservoir by the 
metallic wall of the heat exchanger. The two sections of the refrigeration system are only in thermal contact, but 
mechanically separated. Usually, the heat exchanger is a metallic box filled by the superfluid helium produced in the 
pot. Moreover the thermodynamic properties of the superfluid helium allow to relocate the saturated Hell pot far 
away and to feed the heat exchanger by mechanically decoupled lines (in the TORE Supra system the 1.5 K pot is 
roughly 10 meters apart from the heat exchanger) 92] . 

The pressure fluctuation on the 4.2 K bath are also rejected by the system. The pressure connection between the 
pressurized superfluid reservoir and the saturated helium bath need to be kept as small as possible to improve the 
refrigerator efficiency. 

We remind that no liquid helium flows trough the connections, and that the steady state heat input at the 4.2 K 
HeI-1.7 K Hell interface is produced by thermal conduction in the helium. This heat flow is proportional to the cross 
section of the interface: the smaller the channel, the better is the efficiency of the refrigerator. Usually [5] the pressure 
connection is a tiny tube of about 1 mm cross section or less. Pressure fluctuations in the saturated 4.2 K Hel bath 



are greatly reduced, at the superfluid helium side by the large hydraulic impedance of this channel 7 . 

The last source of noise in the helium bath is the background mechanical noise (seismic and man made) coupled to 
the liquid by the motion of the walls of the refrigerator. By this way the mechanical background noise can spoil the 
rejection of the mechanical-suspension attenuator bypassing the filter via the acoustic wave transmission of the liquid 
helium. To counter this effect we foresee the modification of the refrigerator outlined in fig. the cavity helium 
II reservoir is hydraulically decoupled from the upper part of the pressurized Hell bath using an hydraulic low pass 
filter giving an attenuation of 200 dB at 4 KHz. The needed attenuation can be obtained using a 4 cell low pass filter 
with a cut-off frequency of ~ 4 Hz, producing an asymptotic attenuation of ~ 60 dB per cell at 4 kHz. 

In this way all the mechanical noise, converted to acoustic wave in liquid helium, is confined in the upper reservoir. 
The upper reservoir is also the starting point of the last stage of the mechanical attenuator used to decouple the 
background mechanical noise from the active part of the detector (the resonant RF cavity) and can be (eventually) 
integrated with the mechanical filter. 

The heat transport in the superfluid helium being given by thermal conduction, not by fluid convection, the only 
requirement on the hydraulic filter is to guarantee a mean cross section of the channel large enough to allow a 
temperature drop in the 0.1 K range (worst figure). 



FIG. 18: Modified refrigerator scheme including the hydraulic attenuator needed to decouple the detector from the mechanical 
noise converted to acoustic waves in the cooling bath. 



Usually, (as an example in the refrigerators used for magnet tests) the tiny tube is not even used, the pressure connection being given 
by a non helium tight mechanical separation between the two reservoirs. 
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B. Mechanical dissipation in a liquid helium bath 



The quality factor of a mechanical resonator is given by the ratio of the elastic energy stored in the mode of vibration 
and the power dissipated in one period: 

Q = u-^—, (6.2) 

where u is the resonant frequency. In vacuum, power can be dissipated only by the intrinsic losses of the resonator 
(Pintr) and the resulting intrinsic quality factor is defined as Qo = ujU/Pmtr- For niobium at low temperature and 
few kHz, Qo ~ 10 8 . 

If the resonator is immersed in a fluid, its oscillatory motion causes a periodic compression and rarefaction of the 
fluid near it, and thus produces sound waves. The energy carried away and dissipated by these waves is supplied from 
the kinetic energy of the body. Thus additional losses will be present and the quality factor will be given by: 

Q = ^ P . Zp — (6 - 3) 

*intr i -L fluid 

In a closed vessel only stationary sound waves can exist. For such waves three main dissipation mechanisms can be 
identified [Bf|: 

1) Power losses in the bulk of the fluid due to the classical processes of viscosity and thermal conduction: 
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Power losses in the viscous layer at the solid-fluid boundary, due to the tangential component of the fluid 
velocity at the boundary. The tangential fluid velocity should be zero at the wall, hence a tangential-velocity 
gradient must occur in the boundary layer of fluid, resulting in a viscous dissipation of energy: 



1 vise 2 

3) Power losses in the thermal layer at the solid-fluid boundary. In a sound wave, in fact, not only the density 
and the pressure, but also the temperature, undergo periodic oscillations about their mean values. Near a solid 
wall, therefore, there is a periodically fluctuating temperature difference between the fluid and the wall, even if 
the mean fluid temperature is equal to the wall temperature. At the wall itself, however, the temperature of the 
wall and the adjoining fluid must be the same. As a result a temperature gradient is formed in a thin boundary 
layer of fluid. The presence of temperature gradients, however, results in the dissipation of energy by thermal 
conduction: 

pic 2 ( ujk \ Y 12 7 - 1 f 2 

Ptherm = ~Z ^ / « dS (6.6) 

2p \2p c v J 7 J s 

The calculation of the total energy dissipated into the fluid requires the evaluation of the integrals over the fluid 
volume in eq. 16.41 and over the fluid-solid boundary in eq. 16.51 and eq. 16.61 which, in turn, depend on the details 
of the fluid velocity pattern. In the following section the calculation of the sound wave excited in a closed spherical 
vessel by an oscillating spherical shell is done. 



1. Forced fluid motion in a spherical vessel 



The problem of the forced motion of a fluid in a closed (rigid) vessel can be analytically solved only in simple 
geometries. Here we shall focus on the problem of the calculation of the sound wave excited by an oscillating spherical 
shell of radius a in a fluid which surrounds the shell and is enclosed in a spherical vessel of radius R (see. fig. 11911 . 

A sound wave is completely described by a scalar function cf>(t, r, 9, (j>) (the velocity potential) such that the velocity 
field is given by u = V^>. The mathematical problem can be castcd in the following way: let us find the scalar function 
cj)(t, r, 9, cj>) which satisfies the wave equation (in spherical coordinates) 



Symbol Property name 


Pa 


equilibrium density 


Po 


equilibrium pressure 


c 


sound velocity 


u 


fluid velocity 


n 


shear viscosity 


C 


kinematic viscosity 


K 


thermal conductivity 


Cv 


specific heat at constant volume 


c p 


specific heat at constant pressure 


7 


specific heat ratio c p /c v 



TABLE I: Fluid properties and notation 




FIG. 19: Sketch of the spherical shell (external radius a) enclosed in a concentric spherical vessel (radius R). 



and the boundary conditions: 

(V0-n)| Q = uo/(M) (6-8) 
(V0-n)U = (6.9) 

If we seek a solution of the wave equation for the velocity potential which is periodic in time, having the form: 
4>(t, r, 9, 4>) — <po (r, 6*, cjj) exp (— iut), then we find for <fio the equation: 

V 2 0o + k 2 (j) = (6.10) 

where k 2 = (uj/c) 2 . The most general solution of eq. I6.lt)! is given by the function 

oo I 

= E E Ri(r)Yr(6^) (6.11) 
;=o m=— i 

where Y{ n {9,(f)) are the spherical harmonics, and Ri{r) satisfies the Bessel equation 

R'l + 2 -R[ + (k 2 - U{n ^ 1) ) Ri = (6.12) 



Temperature 




T 


1.5 K 


Pressure 




Po 


10 5 Pa 


Density 




Po 


1.47 x 10 2 kg/(m 3 ) 


Shear viscosity 




V 


1.5 x 10" 6 kg/(m-sec) 


Kinematic viscosity* 




c 


1.5 x 10~ 5 kg/(m-sec) 


Thermal conductivity" 


K 


0.1 W/(m-K) 


Specific heat (const. 


vol.) 


C„ 


1150.3 J/kg 


Specific heat (const. 


press.] 


Cp 


1151.1 J/kg 


Specific heat ratio (c 


p/c v ) 


7 


1.000695 


First sound velocity 




c 


242 m/sec 



"HI sec. 25] 

TABLE II: Physical properties of the liquid helium bath 



External radius 


a 


0.102 m 




Thickness 


w 


2 x 10~ 3 


m 


Density 


PNb 


8.6 x 10 3 


kg/(m 3 ) 


Mass 


M 


2.2 kg 




Angular frequency 


Ul 


2.5 x 10 4 


rad/sec 


Intrinsic quality factor 10 s 






Vessel radius 


R 


0.3 m 




Helium volume 


V 


0.1 m 3 





TABLE III: Characteristics of the small-scale system 



which has solutions: 

Ri{r) = -.-, v - , ■ (g 13) 



AJ^l + ±,kr) + B t J 2 {l + ±,kr) 



where J\ and J 2 are the Bessel functions of the first and second kind. 

If we suppose that the vibrating spherical shell oscillates with quadrupolar pattern (I — 2, m — 2), so that the 
boundary condition at r = a can be put in the form (V<f> ■ xi)\ a = uqF 2 2 (0, <fi) we get the following solution for the 
sound wave induced in the fluid: 

<po(r, 0, (j>) = 7= Y 2 (9, <p) 6.14 

with the two constants A 2 and B 2 determined by the boundary conditions. 



2. Numerical estimates 



Let us focus on a specific example. If we take the radius of the vibrating spherical shell a = 0.1 m and the radius 
of the external vessel R = 0.3 m, we get Q ~ 10 6 , for a niobium sphere, oscillating at / = 4 x 10 3 Hz (the frequency 
of the fundamental quadrupolar mode), in a bath of subcooled, superfhiid helium at T — 1.5 K and P = 10 5 Pa. The 
physical properties used in the calculation are summarized in table ITTT1 

For a large system (a = 0.6 m, R = 1 m), oscillating in the fundamental quadrupolar mode at / = 800 Hz in a 
liquid helium bath with the same characteristics, we obtain Q ~ 10 7 . The physical properties used in the calculation 
are summarized in table llVl 

Obviously, in addition to the mechanism of absorption which have already been considered, energy may also be lost 
from the sound field in the following ways: 

1. by transmission of mechanical energy to the material of the vessel wall; 



External radius 


a 


0.6 m 




Thickness 


w 


0.2 m 




Density 


PNb 


8.6 x 10 3 


kg/(m 3 ) 


Mass 


M 


5.5 x 10 3 


kg 


Angular frequency 


Ul 


5 x 10 3 rad/sec 


Intrinsic quality factor 10 8 






Vessel radius 


R 


1 m 




Helium volume 


V 


3.3 m 3 





TABLE IV: Characteristics of the large-scale system 



2. by subsequent radiation of energy from the outer surface of the vessel to the surrounding medium; This loss 
might be expected to acquire significance at resonant frequencies of the vessel; 

3. by propagation of sound through the liquid helium tubing system. 

These items need to be studied in detail on a realistic model of the whole cryogenic system and, at the end, will be 
the subject of experimental check. 

VII. DETECTION ELECTRONICS 

The three main functions of rf control and measurement system of the experiment are: 

1. The first task of the system is to lock the rf frequency of the master oscillator to the resonant frequency of 
the symmetric mode of the cavity and to keep constant the energy stored in the mode. The frequency lock of 
the master oscillator to the cavity mode is necessary to fill in energy in the fundamental mode of the cavity. 
The reduction of fluctuations of the stored energy to less than 100 ppm greatly reduces the possibility of 
ponderomotive effects due to the radiation pressure of the electromagnetic field on the cavity walls and helps to 
minimize the contribution of the mechanical perturbations of the resonator to the output noise. The frequency 
lock allows also to design a detection scheme insensitive to fluctuations of the resonant frequency of the two 
cavities forming our detector. 

2. The second task is to increase the detector's sensitivity by driving the coupled resonators purely in the symmetric 
mode and receiving only the rf power up-converted to the antisymmetric mode by the perturbation of the cavity 
walls. This goal can be obtained by rejecting the signal at the symmetric mode frequency taking advantage of 
the symmetries in the field distribution of the two modes. Our system, despite of some additional complexity, 
guarantees the following improvements over the one used in previous experiments: 

• a better rejection of the phase noise of the master oscillator obtained using the sharp resonance (Q = 10 10 ) 
of the resonator as a filter; 

• a better insulation of the drive and detector ports obtained by using separate drive and detection arms of 
the rf system; 

• the possibility of an independent adjustment of the phase lag in the two arms giving a better magic-tee 
insulation at the operating frequencies; 

• a greater reliability for the frequency amplitude loop using the transmitted power, instead of the reflected, 
coming from the cavity. 

3. The third task is the detection of the up-converted signal achieving the detector sensitivity limit set by the 
contribution of the noise sources j^]. Slow pressure fluctuations on the cooling bath, hydrostatic pressure 
variations due to the changes in the helium level, pressure radiation, and so on, will change in the same way 
the resonant frequency of both modes. Using a fraction of the main oscillator power as local oscillator for the 
detection mixer, the detection system becomes insensitive to frequency drifts of both modes, allowing for a 
narrow band detection of the up-converted signal produced by the cavity wall modulation. 



A. The rf control loop 



The RF signal generated by the master oscillator (HP4422B) is fed into the cavity through a TWT amplifier giving 
a saturated output of 20 Watt in the frequency range 2-4 GHz. The energy stored in the cavity is adjusted at the 
operating level by controlling the output of the master oscillator via the built-in variable attenuator. 

The output signal is divided by a 3 dB power splitter. The A output of the splitter is sent to the TWT amplifier, 
the B output is sent, through the phase shifter (PS), to the local oscillator (LO) input of a rf mixer acting as a phase 
detector (PD). The output of the rf power amplifier is fed to the resonant cavity through a 

double directional coupler, and a 180° hybrid ring acting as a magic-tee. The rf power enters the magic-tee via 
the sum arm, E, and is split in two signals of same amplitude and zero relative phase, coming out the tee co-linear 
arms 1 and 2. 

The rf signal, reflected by the input ports of the cavity, enters the magic-tee through the co-linear arms. The two 
signals are added at the S arm and sampled by the directional coupler to give information about the energy stored in 
cavity allowing for the measurement of the coupling factor, quality factor, stored energy. While driving the cavity on 
the symmetric mode no reflected signal is shown at the A port of the magic-tee where the signals coming from the 
co-linear arms are algebraically added to zero due to the 180° phase shift. 

To get the maximum of the performances of the magic-tee we need to have equal reflected signals (phase and 
amplitude) at the cavity input ports. To preserve the signal integrity we use matched lines (in phase and amplitude) 
inside the cryostat. Because the phase shift is very sensitive to temperature inhomogeneities between the two cables 
and the phase difference between the two co-linear arms of the magic-tee gives a quite strong signal at the A port, we 
need to compensate for differential thermal contractions of the cables inside the cryostat, leading to phase unbalance 
in the feed lines. To do that we insert a phase shifter in one of the lines to compensate for differences and to reduce 
to a minimum the leakage of the unwanted modes on the two ports. As we will show in section TVII Bl mode leakage 
of the symmetric mode to the A port sets a limit to the system sensitivity increasing the overall noise level of the 
detector. 

Mode leakage of the antisymmetric mode to the E port reduces the system sensitivity by reducing the signal level 
available for detection. The output ports of the cavity are coupled for a maximum output signal on the antisymmetric 
mode (detection mode) and the magic-tee is used to reject the rf power at the frequency of the symmetric mode. A 
fraction of the signal at the E port is fed to the rf input of the phase detector PD via a low noise rf amplifier. The 
intermediate frequency (IF) output of the phase detector PD is fed back to the rf master oscillator to lock the output 
signal to the resonant frequency of the resonator. The total phase shift around the loop is set through the phase 
shifter PS, to have the maximum energy stored in the detector. A careful design of the servo loop amplifier (SLA) 
guarantees the stability of the system and the rejection of the residual noise of the master oscillator up to one MHz. 
The same fraction of the E output of the output magic-tee is used to keep constant, to 100 ppm, the energy stored in 
the cavity feeding back an error signal to drive the electronically controlled output attenuator of the master oscillator. 

Great deal of care is needed in tailoring the frequency response of both controls because the two loops can interact 
producing phase-amplitude oscillations in the rf fields stored in the cavities. 

B. Sensitivity enhancement using the mode symmetry 

The two modes of the detector cavity have (as in the case of two coupled pendula) opposite symmetries of the fields. 

Using two separate sets of ports to drive the cavity and to receive the up-converted signal at the frequency of the 
antisymmetric mode, the cavity acts as a very sharp filter (due to the high Q), with an high rejection of the noise 
coming from the master oscillator at the frequency of the up-converted signal. This already low residual noise, can 
be even more reduced in our scheme using two magic-tees to drive the cavity purely in the symmetric mode and to 
detect only the up-converted energy, rejecting the unwanted field components by an amount given by the magic-tee 
insulation. 

In the case of an ideal magic-tee the mode rejection is infinite. If the cavity is driven purely in the symmetric mode 
no symmetric mode component is transmitted through the system and there will be no signal at the outp ut port. In 
the ideal case this result is obtained also in the more simple scheme used by Melissinos and Reece [22ll3> measuring 
the up-converted power coming out of the detector along the input lines. 

Our scheme gives better sensitivity and performances in the real case. The first obvious gain is the sum of the 
A and E port insulation of the two tees, plus the possibility of adjusting separately the input and output lines to 
get better mode rejection. In a commercial magic-tee the insulation is specified to be « 25 - 35 dB over its own 
bandwidth. The reason for this quite low insulation is mainly due to the difficulty of balancing on a large range of 
frequency the phases of the signals coming from the two co-linear arms of the tee. 

A phase unbalance as small as five degrees reduces the insulation from A to E port to only 25 dB. 



Our electronic scheme allows for an independent compensation of the magic-tee phase mismatch both at the feed 
frequency and at the detection frequency in a flexible way: the phase mismatch is compensated using a variable phase 
shifter at the input of one of the co-linear arms. The optimum phase at the input side results in a pure excitation of 
the symmetric (drive) mode, keeping the power at the frequency of the antisymmetric (detection) mode 70 dB below 
the level of the drive mode. Adjusting the phase at the output will couple only the antisymmetric mode component, 
rejecting the symmetric mode component by 70 dB. The total symmetric-antisymmetric mode rejection of the system 
is the sum of the attenuation we can obtain from the two 180° hybrids. 

The input and output ports of the two cell cavity need to be critically coupled (/3 = 1) to the rf source and to 
the rf detection system. In this way we have the optimum transfer of power to the symmetric mode (a maximum of 
stored energy) and to the antisymmetric mode (a maximum in the detector output). Because the frequency and field 
distribution of the two modes are quite close, the input and output ports are critically coupled to both modes. For 
that reason 50% of the symmetric mode signal is coupled to the idle S port at the output magic-tee, 

and symmetrically 50% of the antisymmetric mode signal is coupled to the idle A port at the input magic-tee. We 
remark that, since those two ports are not used in our detection scheme, it is necessary to reflect back into the cavity 
the energy flowing out of them. This task is achieved by a proper termination of the ports; a careful analysis showed 
that closing the ports with an open circuit completely decouples the input and output arms (optimum rejection of 
the symmetric mode) and maximizes the stored energy and the detector sensitivity. 

At the £ port of the detection arm we insert a directional coupler to sample a tiny amount of the symmetric mode 
power coming from the cavity. This signal is fed into the frequency-amplitude servo loop used to lock the master 
oscillator to the cavity frequency and to keep constant the energy stored in the cavity. 

C. Detection of the converted signal 

Having devised a means of suppressing the high level pump frequency signal, the remaining problem is to amplify 
and detect the low level parametric mode signal. The detection limit, hence overall sensitivity of the machine is now 
determined by the RF thermal noise level. In order to detect the low level signal it must be amplified and to do this 
without degrading the signal to noise ratio, a low noise amplifier is required. The lower the effective noise figure of 
the RF amplifier, the higher will be the sensitivity. In this case we can make use of the fact that the resonant cavities 
are operated at cryogenic temperatures and use a cryogenic low noise amplifier with lowest possible noise figure. 

Since the late 1970s research has been conducted into the performance of GaAs FET and later GaAs HEMT RF 
amplifiers at cryogenic temperatures. The principle applications for this technology have been for radio astronomy 
receivers and for satellite ground station receivers. A typical receiver uses a 2 stage closed cycle Helium refrigerator 
reach a base temperature of around 15 K. HEMT amplifiers work well at this temperature if carefully design, but can 
also work down to liquid Helium temperatures. However, the additional benefit in noise temperature performance of 
operating at LHe compared to 15 K is not great and is mostly not justified when traded-off against of system complex- 
ity, when the equipment must operated in a radio astronomy telescope or a satellite ground station antenna. However, 
in this application, cooling to 4.5 K or 1.5K in liquid Helium is available in the same dewar as the superconducting 
coupled cavities. 

The layout of the detection electronics is shown in fig. 1201 
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FIG. 20: Electric circuit of the detector system 

The signal converted to the antisymmetric mode by the interaction between the mechanical perturbation and the rf 
fields is coupled to the A port of the detection arm of the rf system and amplified by the low noise rf amplifier LNA. 

The converted rf signal amplified by the LNA is fed to the rf input port (RF) of a low noise double balanced mixer 
Ml; the local oscillator port (LO) of Ml is driven by the symmetric mode rf power (at angular frequency u> s ). The 



LO input level is adjusted to minimize the noise contribution from the mixer. As shown in the previous section the 
input spectrum to the RF port of the mixer consists of two signals: the first at angular frequency lo s coming from 
the rf leakage of the symmetric mode through the detection system; the second is the converted energy at angular 
frequency u> a . Both signals are down-converted by the Ml mixer giving to the IF port a DC signal proportional to 
the symmetric mode leakage and a signal at angular frequency f2 proportional to the antisymmetric mode excitation. 
The down-converted IF output is further amplified using a low noise audio preamplifier. 

For the detection electronics of mechanically coupled interactions, at angular frequency f2, as gravitational waves, 
since the exact frequency and phase of the driving source is not known, we can't perform a synchronous detection. 
We need to perform an auto-correlation of the detector output, 

or to cross-correlate the outputs of two different detectors. 

The outlined detection scheme gives the benefit of being insensitive to perturbations affecting in the same way the 
frequency of the two modes. 

The first point to note in terms of the general architecture of the solution is that if a cryogenic amplifier is to 
be used within the same cryogenic dewar as the coupled cavities, then the phase shifter, variable attenuator and 
180° hybrid also need to be inside the dewar. If this is not the case then the noise temperature contribution of 
these components will dominate and there will be little, if any, benefit in having a cooled amplifier. As it is known 
that HEMT amplifiers can be produced for cryogenic operation, we shall first focus on the passive components. A 
consequence of placing the phase shifter and variable attenuator inside the cryogenic dewar is that they will have to 
be controllable by an electrical signal, rather than manually. For this reason the description below considers voltage 
controlled devices. In the interest of implementing a solution which is reasonably economic, the preferred approach 
will be to use commercial devices for the passive components, rather than developing special-to-type components. 
Many electronic components, particularly passive ones, will operate at cryogenic temperatures, even if not specified to 
do so by their manufacturers. Hence, the following focuses on the potential problems which will need to be investigated 
in order to qualify commercial components. 

The 180° hybrid is essentially a just a strip line device. The only concern about operation in cryogenic temperatures 
is due to differential thermal contraction of the materials of the component, particularly in the area of the connections 
between the strip line and coaxial I/O connectors. There may also be changes in dielectric properties of the substrate, 
which may affect the impedance of the I/O ports. Nevertheless strip line couplers have been used in many cryogenic 
applications without any problem, so we shall only validate this component with a simple test. 

Voltage controlled phase shifter. The key element of the shift shifter is a varactor diode. Silicon and GaAs 
diodes are know to be capable of operation at cryogenic temperatures. However, particularly in the case of silicon, 
the electrical characteristics change significantly at cryogenic temperatures. Therefore, it is predicted that a phase 
shifter will still operate a cryogenic temperatures but its transfer characteristics will change significantly from the 
ambient temperature specifications. For this reason a particular aim will be to calibrate the device when operating 
at cryogenic temperatures. Another factor to note is that varactor diode phase shifters have limitations in operating 
power. Between dBm and + 5dBm input power the varactor diodes begin to work as rectifiers. In any case, the 
pump energy in the cavity output ports is typically below OdBm, but this should be confirmed. 

Voltage controlled variable attenuator. The key element of a voltage variable attenuator is a diode and so the 
analysis given above for the phase shifter is also relevant for the variable attenuator. 

Bypass Switch. The function of the bypass switch is to allow coarse tuning of the circuit for pump frequency 
suppression. Before the phasing is tuned the pump frequency signal may be large and could saturate the low noise 
amplifier. The bypass switch will allow for direct connection of a spectrum analyser outside the dewar for initial 
tuning. The bypass switch will be coaxial relay type. 

VIII. RESPONSE OF THE DETECTOR 
A. Signal 

The interaction between the stored e.m. field and the time-varying boundary conditions depends on both how the 
boundary is deformed by the external perturbation and on the spatial distribution of the fields inside the resonator. 
Any field configuration inside the resonator can be expressed as the superposition of the electromagnetic normal modes 
0: E(r, t) = J2 £n(t) E„(r); H(r, t) = £ H n {t) H„(r), with £ n (t) = E • E„ dV; H n (t) = ^ J H • H„ dV and 

J H„ • H m dV = J E„ ■ E m dV = S nm . 

Similarly if u(r, t) denotes the displacement of the mass element of the walls at point r and time t, relative to its 
position at rest, the displacement of the walls can be written as the superposition of the mechanical normal modes 
£a( r ) j2£g: u(r, t) = €a( r ) Qa(t), with J £ a (r) • £p(r)p(r) dV = M5 a p, where M and p are the mass and the density 
of the walls. 



We want to study the energy transfer between two electromagnetic normal modes of a resonator whose boundary 
is perturbed by an external force when mode " 1" is highly excited by an external source at its eigenfrequency, 
Til ~ A\ cos(wii) and E\ ~ A\ sin(wit), with constant amplitude A\. We shall assume that only one mechanical mode 
couples to the external force. Including empirical damping terms, the equations of motion for the field amplitudes 
are given by: 

H 2 + ^H 2 + u\H 2 = -\<4q m C%Hx (8.1) 



^^ m + < qm = ^-\fn 2 ni (8. 2 ) 

where f m {t) = J f (r, t)-£ m (r) dV, and f (r, t) is the volume force density which acts on the walls. The time-independent 
coupling coefficient is given by (the superscript m labels the mechanical normal mode, while the subscripts label 
the electromagnetic modes): 

C% = [ (H a • H x - E 2 ■ Ei) £ m -dS. (8.3) 
Js 

The dependence of the coupling coefficient, and therefore of the energy transfer, both on the field spatial distribution 
and on the boundary deformation, has been checked using a resonator made up of two pill-box cavities, mounted 
end-to-end and coupled by a small circular aperture in their common endwall. The perturbation of the resonator's 
boundary was induced by two piezoelectric crystals mounted in the center of the two circular endwalls. The TEqh 
symmetric mode at 3 GHz was excited by an external rf source and the piezos were driven by a synthesized oscillator 
tuned at the frequency corresponding to the symmetric-antisymmetric mode separation (w 500 kHz). The relative 
phase of the signals driving the two piezos could be set to degrees and to 180 degrees with an external switch. Eq. 
18.31 predicts for this field and for this boundary configuration = 0, for the in-phase excitation and = 2, for 
the excitation with 180 degrees phase lag. These predictions are clearly confirmed by the data shown in Fig. |^ 
The explicit calculation of the coupling coefficient CJJ for an arbitrary deformation of the resonator volume can be 
done only by numerical methods. Analytic calculations showed that for an ideal spherical hollow resonator, excited 
in the fundamental quadrupolar mechanical mode and in the TEon electromagnetic mode, C%{ = 1.9. More detailed 
calculations, made on a realistic model of the coupled spheres, including the central coupling cell and the e.m. input 
and output ports, were made by finite element methods. 

H2(t) is given by: 

W 2 (t) = / K(f)f m (t-T)dr=— / K(uj)f m (u;)exp(iLot)dw, 



K{t) = for t < . (8.4) 



The density of the external forces driving the motion of the system, i.e. the tides generated by the passing gravitational 
wave, are given by 

f i {r,t)=-p{T)R 0iaj {t)x^ (8.5) 

where p is the wall density, and i?oiOj (t) are the components of the Riemann tensor evaluated in the center of mass 
of the solid. The components of the Riemann tensor can be expressed in terms of the adimensional amplitude of the 
g.w. R 0l0j (t) = -l/2oy(i). 

For a "+" polarized plane g.w travelling along the z axis the force density, in the proper reference frame attached 
to a detector lying in the xy plane, has the form: 

f(r,t) = ±p(r)a(t)[-x,y,0] , (8.6) 

with a similar expression for the "x" polarization. The generalized force, acting on the m mechanical mode, then 
has the form 

f m (t) = -\a{t) [ [(U)xX-(U) y y]p(r)dV. (8.7) 

z JVol 

f m (w) is then given by: 



f m (uj) = - MLU 2 £a(td) 



(8.8) 
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FIG. 21: Power trasfer between the symmetric and the antisymmetric mode with in-phase piezo excitation (solid) and 180 
degress out-of-phase excitation (dashed) 



having denned the detector's effective length (for this mechanical mode and wave polarization): 



C = T-fi ltirn)* X -(Zrn) y y}p(r)dV (8.9) 
m -'Vol 



The function K{uS) is found to be: 



~ s eu, 2 A 1 /(2¥) 

*M = 7—- ~ — - ; J\ 1 i JT^j^i *M^it) (s.id) 

[w| - (w + uJi) 2 + i ^ (u + wi)] 



(C^Aio; 2 )2 
4M 



The second term in the denominator of Eq. 18.101 is due to the back-action of the electromagnetic fields on the 
cavity walls. Its effect is particularly important when LO2 — lo\ ~ u) m . In this case optimum signal transfer does not 
correspond to the maximum field amplitude A\ that can be stored in the initially excited mode (that is limited by 
the critical field of the superconductor, H c « 0.2 T for niobium). Instead, optimum signal transfer is obtained from a 
lower value that has to be optimized according to the prevailing experimental conditions and detector characteristics. 
The dependence of the signal on the energy stored in the initially excited mode will be further discussed in Section 
IVIIBI where the expected sensitivity of the detector in various experimental conditions is analyzed. 

The average energy stored in mode number 2 is U2 = l/2\H.2(t)\ 2 and the power extracted from a cavity port 
coupled to an external load with a coupling coefficient /?2 is P2 = P2 {^2/ 2)^2- 

If a(t) = ha(t) with a(t) periodic at angular frequency tt, the expression of the output power from mode 2 is given 
by P 2 = U(Sl) h 2 , with 

ft(fi) = l^QHl 2 C 2 \K(Q)\ 2 . (8.11) 



B. Noise 



To study how the thermal fluctuations of the walls contribute to the output signal, we will start again from Eqs. 
18.1118.21 taking now the external force f m (t) as a stochastic force with constant power spectrum 5/ m = 4MkE.Tu) m /Q m 
[37| . The output noise power spectral density is given by: 

S PP (w) = ^\K(u,)\*S ff (8.12) 

The thermal fluctuations of the electromagnetic field itself can be calculated in an analogus way. Starting again 
from Eqs. I8.1H8.2I we consider an external, stochastic force coupled to the magnetic field H.2(t), with constant power 
spectrum Sf e = 4kT/(uj2Q). The output noise spectral density is given by: 

S H H = ^\x(")\ 2 Sf e (8.13) 
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with x( u; ) given by: 
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*M = T— — " . r — i 77^7- V, • ( 8 - 14 ) 
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We note that in the limit of vanishing coupling, C\\ ~ 0, the contribution to the output noise of the walls vanishes, 
while the contribution of the electromagnetic field fluctuations tends to be like a simple harmonic oscillator of frequency 
LU2 and quality factor Q. 

Other noise sources must also be taken into account. To operate our device we have to feed microwave power into 
mode 1, using an external rf source locked on mode 1, at frequency uj\. The master oscillator phase noise is filtered 
through the resonator linewidth; the power spectral density has the following frequency dependence |35|: 

4ftP,/K8) 

^mo(w) - " — (8.15) 

( J-\ -L I u+Ul _ "1 

where Pi is the power input level and ft is the coupling coefficient of mode 1 to the output load. From the above 
equation we can estimate the microwave power noise spectral density at the detection frequency u) = 0J2 — <-^i '■ 

Smo(co 2 - wi) ( — ) (8-16) 

This figure can be improved if the receiver discriminates the spatial field distribution of the e.m. field at frequency 
u>2, i.e. if it is sensitive only to the power excited in mode number 2, rejecting all contributions coming from mode 
number 1. In this way mode 1 is decoupled from the output load and [3\ — 0. The experimental set-up, based on the 
use of two magic-tees which accomplishes this issue, is discussed in detail in 0. Of course, the mode discrimination 
cannot be ideal, and some power leaking from mode 1 to the detector's output will be present. Nevertheless our 
previous work has demonstrated that with careful tuning of the detection electronics we can obtain ft < 1CT 14 @. 

The input Johnson noise of the first amplifier in the detection electronics has to be added to the previous contribu- 
tions to establish the overall noise spectral density. It can be described by the frequency independent spectral density 

M 

Sjj = k B T{10^ N l^ - 1) = k B T eq (8.17) 

where TV is the noise figure of the amplifier (in dB) and T the operating temperature. The equivalent temperature 
(or noise temperature) T eq , is equal to the temperature (in Kelvin) of a 50 ohm termination at the input of an ideal 
noiseless amplifier with the same gain and generating the same output noise power. 
Other noise sources (e.g. the seismic noise) are not considered in this paper. 

We shall characterize the noise in our detector by a frequency dependent spectral density Sh{ui), with dimension 
Hz -1 defined as 

*M = H (8-18) 
where S n = Spp + Sum. + Smo + Sjj, is the detector noise spectral density. 



C. Sensitivity 

Let us focus our attention on the system based on two spherical niobium cavities working at u>x w 0J2 ~ 2 GHz with 
a maximum stored energy in the initially excited symmetric mode of U\ « 10 J per cell (corresponding to a maximum 
surface magnetic field H max = 0.1 T, half the critical field of niobium). This is a small-scale system with an effective 
length of 0.1 m and a typical weight of 5 kg. The lowest quadrupolar mechanical mode is at u m ~ 4 kHz. In the 
following, we shall consider an equivalent temperature of the detection electronics T eq = 1 K. 

A possible design of the detector uses both the mechanical resonance of the structure, and the e.m. resonance. Due 
to the tuning system, the detection frequency can be made equal to the mechanical mode frequency uo m ss 0J2 — uji- 
The expected sensitivity of the detector for 4 kHz is shown in figure 1221 for a mechanical quality 

factor Q m — 10 3 (solid line) and Q m = 10 6 (dashed line). Note that, in the two cases, the optimum sensitivity is 
obtained with different values of stored energy. In both cases the stored energy has been optimized for maximum 
detector bandwidth. When the mechanical quality factor is higher (Q m = 10 6 ) the stored energy has to be maintained 
much under the maximum allowed value. 

When u>2 — wi = u)m, the detector sensitivity is limited by the walls thermal motion. In this case, a lower T eq would 
increase the detection bandwidth. 
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Since our detector is based on a double resonant system (the mechanical resonator and the electromagnetic res- 
onator) it can be operated also for frequencies u>2 — 7^ u> m . At frequencies L02 — u>i < 1 kHz the master oscillator 
phase noise will, in general, completely spoil the system sensitivity, while at frequencies L02 — u>i > 10 kHz the 
noise coming from the detection electronics will dominate. The expected sensitivity of the small-scale detector for 
u>2 — U\ = 10 kHz is shown in figure 1251 

In order to increase the expected sensitivity a large-scale system has to be developed. A possible design could 
be based on two spherical cavities working at u>% w w 2 ~ 500 MHz, with u m ~ 1 kHz. This system could have a 
maximum stored energy of U\ « 1200 J per cell, an effective length of 0.4 m and a typical weight of 2300 kg. With a 
reasonable choice of system parameters one could obtain the sensitivity shown in figure YIM for the double-resonance 
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FIG. 23: Calculated system sensitivity (uj m w 4 kHz, w 2 - Wi « 10 kHz, Q = 10 10 , T = 1.8 K, T eq = 1 K, stored energy 
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case (u>2 — u>\ = u> m ). As in the previous (small-scale) case the energy store in the initially excite mode has been 
optimized for maximum bandwidth, and it has to be much less then the maximum allowed. Also in this case lowering 
T eq corresponds to an increase of the detection bandwidth. 

Obviously the large-scale system could also be used at higher frequencies; in this case a good sensitivity can be 
achieved in a narrow detection bandwidth (see Fig. 1251) - 

It is worth noting that the narrow detection bandwidth is not an unavoidable drawback of the system. Actually 
its value is determined by the coupling coefficient of the antisymmetric mode at the A port of the output magic- 
tee, and can be adjusted changing this coupling. This corresponds to changing (lowering) the quality factor of the 
antisymmetric mode, leaving the quality factor of the symmetric mode unaffected. Of course an increased bandwidth 
corresponds to a lower sensitivity, since the latter is proportional to the antisymmetric mode quality factor. The 
possibility to increase the detection bandwidth is also interesting for other possible applications of this detection 
technique which is based on the parametric frequency conversion between two electromagnetic modes in a cavity (for 
example, in connection with recently proposed detectors based on the dual resonator concept |38l l39l l40|). 



IX. DATA ANALYSIS 



Several aspects of MAGO data analysis need to be investigated. 

The data-analysis for both single and multiple MAGOs will be akin of those used for acoustic detectors m 
view of the narrowband MAGO response. 

Data analysis strategies for MAGOs operating above / ~ 4 ■ 10 3 Hz will be strictly non parametric, in view of the 
present lack of GW source/signal models in this frequency range. 

Optimum MAGO tuning schedules for detecting unknown (sub)millisecond pulsars 0]> E3l wm be obviously a 
most relevant issue to be investigated. 

A "xylophone" of MAGOs tuned in the range between 10 3 Hz and 10 4 Hz might be able to both detect and estimate 
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FIG. 24: Calculated system sensitivity (ui m « lo 2 - Oi ~ 1 kHz, Q = 10 10 , Q m = 10 6 , T = 1.8 K, T eg = 1 K, stored energy 
U ~ 1 J per cavity) 



the chirp-mass Ii3 of galactic light 8 BH-MACHO binaries which would produce a distinct signature in the 
xylophone output |46| . The expected (optimistic) event-rate might be non negligible |4l| • 

Once MAGO prototypes will be up and running, a primary task will be to investigate the experimental noise PSD 
and transfer function, by comparison to their theoretical counterparts. 

X. FUTURE DEVELOPMENTS 

The sensitivity of the device currently under test, consisting of a couple of nearly spherical niobium cells, 2 mm 
thick, whose total mass is approximately 4.5 kg, might be improved in several ways. Among these the most promising 
are: 

• an increase of the mechanical quality factor of the structure; 

• an increase of the electromagnetic quality factor of the superconducting cavities; 

• a lower noise temperature of the detection electronics; 

• the increase of the size (and mass) of the detector; 

• the use of several detectors in an array or in a xylophone. 
All these items are discussed in the following. 



The highest GW frequency radiated during binary inspiral is twice the orbital frequency of the last stable circular orbit, and is l44l 
/ ~ 4- 10 3 (M Q /M) [Hz] . 



10 *3 I , , , 1 , , , 1 , , , 1 , , , i , , , 1 

4001.4 4001.6 4001.6 4002 4002.2 4002.4 



Frequency [Hz] 

FIG. 25: Calculated system sensitivity (u m w 1 kHz, uj 2 - Ui « 4 kHz, Q = 10 10 , Q m = 10 6 , T = 1.8 K, T eq = 1 K, stored 
energy U « 1200 J per cavity) 



A. Increasing the mechanical quality factor 

Although the S/N ratio does depend on several parameters, we can state that if the noise power spectral density 
is limited by thermal noise, the minimum Sh is proportional to the square root of the mechanical quality factor, 

The mechanical quality factor depends both on the intrinsic mechanical energy dissipated in the material of the 
resonator and on the energy dissipated through the coupling to the external environment. The intrinsic quality factor 
of niobium at low temperature and few kHz is Q m ~ 10 8 . Since our detector is immersed in a bath of subcooled 
superfluid helium, improvements in Q m can be expected thanks to carefully designed liquid helium bath and suspension 
system 9 . 

The limit imposed by mechanical dissipation in the superfluid helium bath is Q m ~ 10 6 for the present small-scale 
detector (see sec. IVI Bfl . 

This figure can be further improved for a larger detector [M > 10 4 kg). In this case, the expected liquid-helium- 
limited mechanical quality factor should be Q m ~ 10 7 (see sec. IVI Bfl . 

We remark that the above limits are essentially due to the will of using a well established refrigeration scheme for 
the superconducting cavity. In principle different schemes could be developed which do not make use of a bath of 
coolant in contact with the cavity (e.g. pulse tube refrigerators). Obviously this requires a consistent R&D effort 
which is not, for the time being, in the line of our proposal which is entirely based on already existing technologies. 



In this document we consider this cryogenic scheme. Different schemes, not based on the use of a liquid helium bath (e.g. pulse tube 
refrigerator), might be foreseen. 



B. Increasing the electrical quality factor 



The gain in sensitivity due to the electrical quality factor, Q, is less trivial to discuss. First of all a distinction 
should be made between the electromagnetic quality factor of the symmetric (pump) mode and of the antisymmetric 
(detection) mode 10 . The quality factor of the pump mode, Qi, strongly affects the carrier (oscillator) phase noise 
and the energy stored in the cavity for a given input power. The quality factor of the detection mode, Q2, affects the 
shot noise power spectral density and the detection bandwidth. 

For the sake of simplicity in our proposal these two quality factor were taken as equal: Qi = Q2 — Q- This 
is certainly a reasonable assumption as far as the intrinsic quality factor of the niobium spheres: Qo — G/R s is 
considered. In fact Q depends only on the geometrical shape of the resonator, on the electromagnetic resonant mode 
and on the surface quality. The overall (loaded) quality factor of a resonator, though, is given by Q = Qo/(l + 0), 
where (3 is the coupling coefficient to the external load 11 . 

As a consequence, it is in general convenient to keep the quality factor of the pump mode Qi as high as possible, 
in order to have the minimum power dissipation for a given energy stored in the detector (or the maximum stored 
energy for a given level of power dissipation) and the minimum carrier phase noise at the detection frequency. This 
requirement is fulfilled when j3\ ~ 1 (critical coupling). 

On the other hand, the value of the detection mode quality factor, Q2, is a parameter that has to be set to the 
optimum value (which depends on the specific experimental situation) by tuning the output coupling coefficient /3% . To 
be more specific, since shot noise and amplifier back-action noise are roughly proportional to (Q2) -1 ^ 2 , when either 
of these is the dominant noise contribution an higher Q2 gives better sensitivity (and, in general, a small detection 
bandwidth). On the contrary, when sensitivity is limited by other noise sources (e.g. thermal noise) it is convenient 
to lower Q2 (increasing fa) to open the bandwidth. 

Improvements in the intrinsic quality factor, Qo, can be expected thanks to specific Nb surface treatments that 
give a lower R s . Quality factor values Qo ^ 10 11 have already been measured on several prototypes of elliptical 
accelerating cavities (R s ~ 1 nOhm) both bulk niobium 0], and sputtered niobium on copper [43. Since those 
cavities were run on the TM io accelerating mode, we can expect Q ~ 10 12 (on the TE n mode) with a comparable 
surface quality, due to the increased geometrical factor 12 . 

C. Increasing the performance of the detection electronics 

A further headway can be foreseen using a low noise amplifier. The noise temperature of a cryogenic, commercially 
available, high-frequency amplifier (GaAs or InP HEMT device) is T eq = 2 K. Bearing in mind that we want to amplify 
an e.m. field at a frequency around u)/ (Itt) ~ 2 GHz (much greater than the detection (g.w.) frequency VL gw = LJ2—W1) 
one can easily see that the minimum (quantum limited) noise temperature of an amplifier is T eq — huj/ks ~ 0.1 K. 

This apparently high noise temperature though is not a drawback, since the parametric transducer operates the 
up-frequency conversion from low energy quanta h£l gw — > htv with no added noise. In this process the gain (in energy) 
is cj/n gw ~ 10 6 [H sec. 7.9-7.11]. 

In other words, the noise temperature seen at the input of the parametric up-converter is T eq ~ 0.1 x Q gw /u> ~ 
10~ 7 K. 



D. Increasing the size of the detector 

Starting from the experience gathered with the existing small-scale prototype, a plausible detector (inner radius of 
the spherical cell a = 0.2 m, thickness w = 5 mm) could be designed and built 13 . 



We remind that our device is based on the principle of parametric conversion of power between the two electromagnetic modes of 
the coupled resonators: the symmetric mode at frequency ld\ and the antisymmetric mode at frequency oJ2, induced by an external 
perturbation of the cavity boundary at frequency Q ~ \u>2 — wi 

The coupling coefficient /3 determines the impedance match between the input (for the pump mode) or output (for the detection mode) 
impedance of the cavity Z c (uj), and the impedance of the external load Zq. The ratio Z c /Zq is simply the coupling coefficient j3. 
The ratio between the geometric factor of a spherical cavity operating in the TE011 mode and an elliptical cavity operating in the TM010 
mode is Gte/Gtm ~ 3.4 

There are no obstacles, in principle, to envisage an even larger detector, facing, obviously, major technical challenges. Major differences 
lay in the greater mass, greater energy stored, lower detection frequencies and possibly fabrication methods (e.g. sputtered Nb on a 
substrate). 



Such detector might have the following parameters: overall mass M ~ 44 kg (approx. a ten-fold increase with 
respect to the present small-scale prototype), characteristic length C ~ 1 m, working rf frequencies uj / (2tt) ~ 1 GHz, 
stored energy [7 ~ 210 J, detection frequency range 2 kHz < Sl gt0 /(27r) < 10 kHz. 

Being the detector tunable in frequency, we shall divide its use into two different working ranges: a Near Mechanical 
Resonance range (hereafter NMR) and an Off Mechanical Resonance range (OMR). NMR operation can only be used 
next to both quadrupolar mechanical resonances of the spherical shell (at frequencies /i ~ 2 kHz and fi ~ 8 kHz) . 
The two ranges require quite different parameter sets and lead to optimized use of the detector, depending on the 
sources under investigation. 

Fig. 1261 and fig. 1271 show system sensitivity calculated with Q m = 10 6 (limited by mechanical dissipation in liquid 
helium), Qq = 8.5 x 10 11 and a nearly quantum limited electronics T eq = 0.05 K, in NMR and OMR mode. All other 
detector parameters are optimized for NMR and OMR operation. As expected, the use of T eq = 0.05 K for the noise 
temperature of the amplifier does not bring a dramatic increase in the peak performance of the system, whereas the 
bandwidth is strongly enhanced. 

We remind that the bandwidth of the detector can be adjusted also by means of the electromagnetic coupling (/3a) 
as discussed in the previous section. Adjusting the electromagnetic coupling is of paramount importance because 
it allows trading maximum sensitivity versus bandwidth and is particularly useful when the maximum sensitivity is 
limited by thermal noise (NMR operation). 
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FIG. 26: Near mechanical resonance PSD of one MAGO. Minimum {S h ) 1/2 ~ 6 x 1(T 21 (Hz)~ 1/2 , bandwidth A/ ~ 350 Hz. 
Mechanical quality factor Q m = 10 6 (limited by losses in liquid helium), electrical intrinsic quality factor Qo = 8.5 x 10 11 
(R s = 1 nohm), amplifier equivalent temperature T eq = 0.05 K (quantum limit). 



E. Increasing the number of the detectors 

A possibly nicer extension for the detector discussed in the previous section is given by multi-detector operation 
to achieve low false-alarm probabilities in coincidence operation. Furthermore, since MAGO can be designed to work 
at any chosen frequency in the range 10 3 -10 4 Hz with the same sensitivity, a "xylophone" of MAGOs might be able 
to both detect and estimate the chirp mass of galactic light BH-MACHO binaries, which would produce a distinct 
signature in the xylophone output. 
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FIG. 27: Off mechanical resonance PSD of one MAGO. Minimum {S h ) 1/2 ~ 1.6 x 10" 21 (Hz) _1/2 , bandwidth A/ ~ 1 Hz. 
Other parameters as in fig. 1261 



The expected (optimistic) event-rate might be not negligible (see section IlII CI) . 

Several authors have already pointed out the usefulness of array detectors, both for source localization and for 
improving the overall signal information (multimode transducer) [BlT . 153 . IBl^ . 

The SNR of N optimally oriented detectors, sharing the same sensitivity band is approx. N times the SNR of the 
single detector [El|. A set of 10 MAGOs would therefore have a peak (S h ) 1/2 ~ 1CT 22 (Hz)" 1 / 2 . 

In building a MAGO coincidence network, several features of the small single detector add to the many known 
benefits of arrays: 

a) The expected cost and requirements confer a non negligible freedom in choosing the installation sites and 
financing; 

b) There is a greater freedom in the choice of the material; 

c) Short thermal cycle allow faster development phase and more flexible management of the instrument; 

d) High frequency operation is preferable from an instrumental point of view because many problems like vibration 
isolation, transducer optimization and amplifier matching are easier to solve; 

e) MAGO size makes it possible to more easily achieve optimal orientation among distant detectors. It makes also 
possible to host a local array in a reasonably limited area. This feature is particularly important in view of the 
detection of a stochastic background of GWs (see sec. IIII E() . 

f ) Thanks to the frequency and bandwidth tuning ability, one could have a nearly complete overlap of the sensitivity, 
allowing narrow-band/high-sensitivity detection mode (OMR operation) not be a strong limitation. This feature 
might play a role in the detection of GW signals from binary mergers (see sec. IIII All . 

g) High frequency signal sensitivity is apt to constructing high resolution images (see fig. I28fl . 



FIG. 28: Reconstructed image (simulated) from a network of 10 MAGO detectors, with SNR — 0.5 for the single detector. 
The angular distance between the two sources is A8 — A(/> = 0.1° and the shown portion of the sky is 1 ° square. 
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